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Abstract Pulsars steadily dissipate their rotational energy via relativistic winds. 
Confinement of these outflows generates luminous pulsar wind nebulae, seen across the 
electromagnetic spectrum in synchrotron and inverse Compton emission, and in optical 
emission lines when they shock the surrounding medium. These sources act as impor- 
tant probes of relativistic shocks, particle acceleration and of interstellar gas. We review 
the many recent advances in the study of pulsar wind nebulae, with particular focus on 
the evolutionary stages through which these objects progress as they expand into their 
surroundings, and on morphological structures within these nebulae which directly trace 
the physical processes of particle acceleration and outflow. We conclude by considering 
some exciting new probes of pulsar wind nebulae, including the study of TeV gamma-ray 
emission from these sources, and observations of pulsar winds in close binary systems. 
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1 INTRODUCTION 



Neutron star: a 

compact degenerate 
stellar remnant, 
formed in the 
core-collapse of a 
massive star. 
Pulsar: A rapidly 
rotating. 

highly-magnetized 
neutron star, which 
generates coherent 
beams of radiation 
along its magnetic 
poles. 

Pulsar wind 
nebula: A bubble 
of shocked 
relativistic particles, 
produced when a 
pulsar's relativistic 
wind interacts with 
its environment. 



The Crab Nebula (Fig.^) is a lmost certainly associated w ith a supernova (SN) ex- 
plosion observed in 1054 CE ((Stephenson &: Green 2002, and references therein) . 
However, this source differs substantially from what is now seen at the sites of 
other recent SNe, in that the Crab Nebula is centrally filled at all wavelengths, 
while sources such as Tycho's and Kepler's supernova remnants (SNRs) show a 
shell morphology. This and other simple observations show that the Crab Nebula 
is anomalous, its energetics dominated by continuous injection of magnetic fields 
and relativistic particles from a central source. 

A 16th magnitude star embedded in the Crab Nebula was long presu med to 
be the stellar remnant and central engine ((Minkowski 1P42l . lPacmi 1967I L This 
was confirmed when 33-ms optical and radio pulsations were detected from thi s 
star in the late 1960s ( Cocke. Disnev fc Tavlor 19691 : Staelin fc Reifenstein 196sh . 
and thes e pulsations were then sh own to be slowing down at a rate of 36 ns 
per day ( Richards fc Cornelia 19 69L The conclusion was quickly reached that 
the Crab Nebula contains a rapidly rotating young neutron star, or "pulsar", 
formed in the SN of 1054 CE. The observed rate of spin down implies that kinetic 
energy is being dissipated at a rate of ~ 5 x 10 38 ergs s , a val ue similar t o 
the inferred rate at which energy is being supplied to the nebula ( Gold 19691 ). 
Following this discovery, a theoretical understanding was soon developed in which 
the central pulsar generates a magnetized particle wind, whose ultrarelativistic 
electrons and positrons radiate synchrotron emission across the electromagnetic 
spectrum ( Pacini Sz Salvati 19731 . Rees Sz Gunn 1974h . The pulsar has steadily 
released about a third of its total reservoir of ~ 5 x 10 49 ergs of rotational energy 
into its surrounding nebula over the last 950 years. This is in sharp contrast to 
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Figure 1: Images of the Crab Nebula (G184.6-5.8). (a) Radio synchrotron emis- 
sion from the confined wind, with enhancements along filaments, (b) Optical 
synchrotron emission (blue-green) surrounded by emission lines from filaments 
(red), (c) Composite image of radio (red), optical (green) and X-ray emission 
(blue), (d) X-ray synchrotron emission from jets and wind downstream of the 
termination shock, marked by the inner ring. Note the decreasing size of the syn- 
chrotron nebula going from the radio to the X-ray band. Each image is oriented 
with north up and east to the left. The scale is indicated by the 2 arcmin scale 
bar, except for panel (d), where the 20 arcsec scale bar applies. 



shell-like SNRs, in which the dominant energy source is the ~ 10 51 ergs of kinetic 
energy released at the moment of the original SN explosion. 

Observations over the last several decades have identified 40 to 50 further 
sources, in both our own Galaxy a nd in the Magellanic Clouds, with properties 
similar to those of the Crab Nebula ( Green 2o"o3 : Kaspi. Roberts fc Harding 20061 ) 
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SN(e): 

Supernova(e) 
SNR: Supernova 
remnant 

PWN(e): Pulsar 
wind nebula(e) 
ISM: Interstellar 
medium 

MHD: Magnetohy- 
drodynamic 
IC: Inverse 
Compton 



- these sources are known as "pulsar wind nebulae" (PWNe). 1 Sometimes a 
PWN is surrounded by a shell-like SNR, and the system is termed "composite" 
(see Fig. |2J). In other cases, best typified by the Crab itself, no surrounding shell 
is seen. 2 

More recently, an additional category of PWNe has been identified, in which 
pulsars with high space velocities produce nebulae with cometary or bow shock 
morphologies as they move through the interstellar medium (ISM) at supersonic 
speeds. The sample of such sources is currently small, but high spatial resolution 
observations, especially in the X-ray band, are rapidly adding to this group. 

Because PWNe have a well-defined central energy source and are close enough 
to be spatially resolved, they act as a marvelous testing ground for studying 
both relativistic flows and the shocks that result when these winds collide with 
their surroundings. Studies of PWNe, particularly the spectacular images now 
being taken by the Chandra X-ray Observatory, allow us to resolve details of 
the interaction of relativistic flows with their surroundings that may never be 
possible in other classes of source, and can provide the physical foundation for 
understanding a wide range of astrophysical problems. 

We here review current understanding of the structure and evolution of pulsar 
wind nebulae, with an emphasis on the explosion of new data and new ideas that 

have emerged in the last few years. Our focus is p rimarily observational; theoreti- 

cal considerati ons h ave been recently discussed bv lvan der Swaluw. Downes &; Keegan 
Melatos r20ol ) and IChevalier C2005D . 



The outline of this review is as follows: in ^we explain the basic observational 
properties of pulsars and their nebulae; in we review current understanding 
of the evolutionary sequence spanned by the observed population of PWNe; in 
50] we discuss observations of PWNe around young pulsars, which represent the 
most luminous and most intensively studied component of the population; in 
we consider the properties of the bow shocks produced by high velocity pulsars; 
and in ^Hlwe briefly describe other recent and interesting results in this field. 



2 OVERALL PROPERTIES 
2.1 Pulsar Spin Down 

Since a pulsar's rotational energy, E ro t, is the source for most of the emission 
seen from PWNe, we first consider the spin evolution of young neutron stars. 

2.1.1 Spin-Down Luminosity, Age and Magnetic Field An isolated 
pulsar has a spin period, P, and a period derivative with respect to time, P = 
dP/dt, both of which can be determined from observations of the pulsed signal. 
The "spin down luminosity" of the pulsar, E = —dE ro t/dt, is the rate at which 

PWNe are also often referred to as " plerions". However, given this term's obscure origin 
JShakeshaft 19791. 1 Weiler fc Panagia 19781) . we avoid using this terminology here. 

2 The absence of a shell around the Crab Nebula is presumably because it has not yet i nter- 
acted with sufficient surrounding gas llFrail et al. 19951 : ISeward. Gorenstein fc Smith 20051) . 
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rotational kinetic energy is dissipated, and is thus given by the equation: 

E = 4^I^ (1) 

where / is the neutron star's moment of inertia and is usually assumed to have 
the value 10 45 g cm" 2 . Values of E for the observed pulsar population range be- 
tween ~ 5 x 10 38 ergs s~ 1 for the Crab pulsar and PSR J05 37-6910, down to 3 x 



10 28 ergs s _1 for the slowest known pulsar, PSR J2144-3933 ( Manchester et al. 200F 
Typically only pulsars with E ^ 4 x 10 36 ergs s -1 (of which ~ 15 are currently 
known) produce prominent PWNe ( Gotthelf 20041 ) . 



The age and surface magnetic field strength of a neutron star can be inferred 
from P and P, subject to certain assumptions. If a pulsar spins down from 
an initial spin period Po such that Ct = —k£l n (where = 2ir/P and n is the 
"braking index"), then the age of the system is ((Manchester fe Tavlor 1977 ): 

(2) 



(n - 1)P 



p \ n— 1 



where we have assumed A; to be a constant and n ^ 1. The braking index, n, has 
only been confidently measured for four pulsars ( Livingstone. Kaspi Gavriil 20051 . 



and references therein), in each case falling in the range 2 < n < 3. 

If for the rest of the population we assume n = 3 (corresponding to spin 
down via magnetic dipole radiation) and Po *C Pi Equation @ reduces to the 
expression for the "characteristic age" of a pulsar, 

r c = -i. (3) 
2P y ' 

Equation © often overestimates the true age of the sy stem, indicating that Pq 



3VS' 

is not much smaller than P (e.g., iMigliazzo et al. 20021 ). PWNe resembling the 



Crab Nebula tend to be observed only for pulsars younger than about 20 000 years 
(see fJH; older pulsars with high space velocities can power bow-shock PWNe (see 
©• ' 

In the case of a dipole magnetic field, we find k = 2M|/37c 3 , where M± is the 
component of the magnetic dipole moment orthogonal to the rotation axis. We 
can thus calculate an equatorial surface magnetic field strength: 

B p = 3.2 x 10 19 (PP) 1/2 G, (4) 

where P is in seconds. Magnetic field strengths inferred from Equation © range 
between 10 8 G for recycled (or "millisecond") pulsars up to > 10 15 G for "mag- 
netars" . Most pulsars with prominent PWNe have inferred magnetic fields in the 
range 1 x 10 12 to 5 x 10 13 G. 

2.1.2 Time Evolution of E and P A pulsar begins its life with an initial 
spin-down luminosity, Eq. If n is constant, its spin-down luminosity then evolves 
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with time, t, as fe.g.. lPacini fc Salvati 1973T ): 



(ra+l) 

£ = £o(^l + i-J (5) 

where 

_Po 2r c 

Tn = — = t (6) 

(n-l)Po ""I 

is the initial spin-down time scale of the pulsar. The pulsar thus has roughly 
constant energy output until a time tq, beyond which E oc t~ ("■+ 1 )/( ri ~ 1 ). The 
spin period evolves similarly: 

i 

t \ i-i 



? = *>(! + - , (7) 
so that P « P for i < r , but at later times P oc t 1 ^™ -1 ). 
2.2 Radio and X-ray Emission from PWNe 

As discussed in H4.1( the resultant deposition of energy with time generates a pop- 
ulation of energetic electrons and positrons, which in turn powers a synchrotron- 
emitting nebula. Radio synchrotron emission is characterized by a power-law 
distribution of flux, such that S u oc v a , where S u is the observed flux density 
at frequency v, and a is the source's "spectral index." At X-ray energies, the 
emission is often described as a power-law distribution of photons, such that 
Ne oc E~ t ', where Ne is the number of photons emitted between energies E and 
E + dE, and T = 1 — a is the "photon index." Typical indices for PWNe are 
—0.3 <^ a <^ in the radio band, and T ~ 2 in the X-ray band. This steepening of 
the spectrum implies one or more spectral breaks between these two wavebands, 
as discussed in H4.6I 

If the distance to a PWN is known, the radio and X-ray luminosities, Lr and 
Lx-, respectively, can be calculated over appropriate wavelength ranges. Typical 
ranges are 100 MHz to 100 GHz for Lr and 0.5-10 keV for Lx- Observed values 
for Lr and Lx span many orders of magnitude, but representative values might 
be Lr ~ 10 34 ergs s _1 and Lx ~ 10 35 ergs s _1 . The efficiency of conversion of 
spin-down luminosity into synchrotron emission is defined by efficiency factors 
7?b = Lr/E and rjx = Ly/E. Typical values are vr ^ 10 -4 and r\x ~ 10~ 3 
([Becker k, Triimper 1997 . Frail fc Scharringhausen 19971 ). although wide excur- 



sions from this are observed. Note that if the synchrotron lifetime of emitting 
particles is a significant fraction of the PWN age (as is almost always the case 
at radio wavelengths, and sometimes also in X-rays), then the PWN emission 
represents an integrated history of the pulsar's spin down, and t\r and r\x are 
not true instantaneous efficiency factors. 
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3 PULSAR WIND NEBULA EVOLUTION 

We now consider the phases of evolution which govern the overall observational 
properties of PWNe. The detaile d theoretical underpinning fo r these evolution- 
ary p h ases is given in studies bv iRevnolds fc Chevalier Chevalier 
20051 ) . iBlondin. Chevalier fc Frierson C200lf ). iBucciantini et al. (200$ ) and van 



der Swaluw et al. (2004). 

We defer a discussion of the details of how the wind is generated to £0J and here 
simply assume that the pulsar's continuous energy injection ultimately results in 
an outflowing wind which generates synchrotron emission. 

3.1 Expansion into Unshocked Ejecta 

Since a pulsar is formed in a SN explosion, the star and its PWN are initially 
surrounded by an expanding SNR. The SNR blast wave at first moves outward 
freely at a speed > (5 — 10) x 10 3 km s , while asymmetry in the SN explosion 
gives the pulsar a random space velocity of typical magnitude 400-500 km s~ 1 . 
At early times the pulsar is thus located near the SNR's center. 

The pulsar is embedded in slowly moving unshocked ejecta from the explosion 
and, since t <C tq, has constant energy output so that E « Eq (see Eqn. 0). 
The pulsar wind is highly over-pressured with respect to its environment, and the 
PWN thus expands rapidly, moving supersonically and driving a shock into th e 



ejecta. In the spherically symmetric case, the PWN evolves as (jChevalier 19771 ): 
R PWN « l.BE^E^M-^, 



1.1 pc 



E \ 1/5 f E SN \^°( M ej \'^f t 



10 38 ergs s _1 / \10 51 ergs/ \10 M Q / \10 3 years 



where Rpwn is the radius of the PWN's forward shock at time t, and E$n and 
M e j are the kinetic energy and ejected mass, respectively, released in the SN. 

Since the PWN expansion velocity is steadily increasing, and the sound speed 
in the relativistic fluid in the nebular interior is c/\^3, the PWN remains cen- 
tered on the pulsar. Observationally, we thus expect to see a rapidly expanding 
SNR, with a reasonably symmetric PWN near its center, and a young pulsar 
near the center of the PWN. A good example of system at this stage of evolu- 
tion is the recently discovered pulsar J1833-1034, which powers a bright X-ray 
and radio PWN, which in turn lies at the center of the young SNR G21.5-0. 9 
(Fig. Ha); ICamilo et al. 20ol. Idiota et al. 20051. iMatheson Safi-Harb 2005h. 



This system is estimated to be ~ 1000 years old. 

3.2 Interaction with the SNR Reverse Shock 



As the expanding SNR sweeps up significant mass from the ISM or circumstellar 
medium, it begins to evolve into the "Sedov- Taylor" phase, in which the total 
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Supernova 
Blast Wave 
and Swept-up 
Shell 




Figure 2: (a) A deep Chand ra X-ray image of the composite SNR G21.5-0.9 
( Matheson fc Safi-Harb 200^ . A circular SNR of diameter « 5' surrounds a 



symme tric PWN of diameter 1'.5, with th e young pulsar J1833-1034 at the 



center (|Camilo et al. 200.1 iGupta et al. 200^ ). The central location of the pulsar 
and PWN and the symmetric appearance of the PWN and SNR both argue for a 
relatively unevolved system, in which the PWN expands freely and symmetrically 
into the unshocked interior of the SNR. (b) A schematic diagram of a composite 
SNR showing the swept-up ISM shell, hot and cold ejecta separated by the reverse 
shock, and the central pulsar and its nebula. The expanded PWN view shows 
the wind termination shock. Note that this diagram does not correspond directly 
to G21.5-0.9, in that a significant reverse shock has probably yet to form in this 
young SNR. 



energy is cons erved and is partitioned e qually between kinetic and thermal con- 
tributions (see Truelove fe McKee 19991 . for a detailed discussion). 



The region of interaction between the SNR and its surroundings now takes on a 
more complex structure, consisting of a forward shock where ambient gas is com- 
pressed and heated, and a reverse shock where ejecta are decelerated. The two 
shocks are separated by a contact discontinuity at which instabilities can form. 
The reverse shock at first expands outward behind the forward shock, but eventu- 
ally moves inward. In the absence of a central pulsar or PWN, and assuming that 
the SNR is expanding into a constant density medium (which, gi ven the effects o f 
progenitor mass loss by stellar winds, may not be the case: see IChevalier 2005 ). 



the reverse shock reaches the SNR center in a time (|Revnolds fc Chevalier 1984 ) : 



M e3 E SN \- l ' 2 ( n 



^" 7 'l0MeJ [W^J [l^=* ' kyr ' (9) 



where no is the number density of ambient gas. At this point the SNR interior is 
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entirely filled with shock-heated ejecta, and the SNR is in a fully self-s imilar stat e 
which can be completely described by a small set of simple equations ( Cox 1973 ) . 
The radius of the shell's forward shock now evolves as Rsnr oc t 2 / 5 . 

In the presence of a young pulsar, the inward moving SNR reverse shock col- 
lides with the outward moving PWN forward shock after a time t co u < ts e dov, 
typically a few thousand years (van der Swaluw et al. 2001, Blondin et al. 2001). 
Even in the simplest case of a stationary pulsar, an isotropic wind and a spherical 
SNR, the evolution is complex. The reverse shock compresses the PWN by a large 
factor, which responds with an increase in pressure and a sudden expansion. The 
system reverberates several times, resulting in oscillation of the nebula on a time 
scale of a few thousand years, and a sudden increase in th e nebular magnetic field 
which serves to burn off the highest energy electrons ((Bucciantini et al. 20031 . 



Reynolds fc Chevalier 19841 . van der Swaluw et al. 200l] ). The crushing of the 



PWN produces Rayleigh- Taylor instabilities, which can produce a chaotic, fila- 
mentary structure and mixing of thermal and non-thermal material within the 
PWN (Chevalier 1998, Blondin et al. 2001). 

In a more realistic situation, the pulsar's motion carries it away from the SNR's 
center by the time the reverse shock collides with the PWN. Furthermore, if the 
SNR has expanded asymmetrically, then the reverse shock moves inward faster 
on some sides than on others. This results in a complicated three-dimensional 
interaction, spread over a significant time interval, during and after which the 
PWN can take on a highly distorted morphology and be significantly displaced 
from the pulsar position (Chevalier 1998, van der Swaluw et al. 2004). An 
example of such a system is the Vela SNR, shown in Figure Ofa). 

3.3 A PWN inside a Sedov SNR 

Once the reverberations between the PWN and the SNR reverse shock have faded, 
the pulsar can again power a steadily expanding bubble. However, the PWN now 
expands into hot, shocked, ejecta at subsonic speeds. In the spherically symmetric 
case, there are two solutions, depending on whether t < To or t > tq (see Eqn. 0). 
In the former ca se, E is approximately const ant and the PWN radius evolves as 
Rpwn oc t 11 / 15 (|van der Swaluw et al. 2001 ). whil e for the latter situation E i s 



decaying, and (for n = 3) we expect Rpwn oc i 3 / 10 ( Reynolds fc Chevalier 1984h . 

At this point, the distance traveled by the pulsar from the explosion site can 
become comparable to or even larger than the radius of an equivalent spherical 
PWN around a stationary pulsar. The pulsar thus escapes from its original wind 
bubble, leaving behind a "relic PWN", and generating a new, smaller PWN 
around its current position (van der Swaluw et al. 2004). Observationally, this 
appears as a central, possibly distorted radio PWN, showing little corresponding 
X-ray emission. The pulsar is to one side of or outside this region, with a bridge of 
radio and X-ray emission linking it to the main body of the nebula. An example 
is the PWN in the SNR G327. 1-1.1, shown in Figure E2» . 

The sound speed in the shocked ejecta drops as the pulsar moves from the 
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center to the edge of the SNR. Eventually the pulsar's space mot ion becomes su- 
personic, and it now drives a bow shock throu gh the SNR interior ( Chevalier 19981 : 
van der Swaluw. Achterberg fc Gallant 19981 ). The ram pressure resulting from 



the pulsar's motion tightly confines the PWN, so that the nebula's extent is 
small, <J 1 pc. Furthermore, the pulsar wind is in pressure equilibrium with its 
surroundings, so that the PWN no longer expands steadily with time. 

For a SNR in the Sedov phase, the transition to a bow shock takes place 
when the pulsar has moved 68% of the distance between the center and the 
forward shock of the SNR (van der Swaluw et al. 1998, 2003). The pulsar is now 
surrounded by a Mach cone, and the PWN takes on a cometary appearance at 
X-ray and radio wavelengths. An example of such a system is PSR B1853+01 in 
the SNR W44, as shown in Figure Ufa). 

A pulsar will typically cross its SNR shell after ~ 40 000 years (see Eqn. 
in Sj2J). If the SNR is still in the Sedov phase, the bow shock has a Mach number 




(b) 




RA (J2Q00) 



15 h 55™30 1 55 ra 00 B 54 m a0 5 54 m 00' 53™30'' 
Right Ascension U2000) 



Figure 3: fa) A 2.4-GHz Parkes map of the Vela SNR (G263.9-3.3), 
( Duncan et al. 199fil ). A limb-brightened shell and a central radio PWN can both 
be seen. The cross indicates the location of the ass ociated pulsar B083 3-45, while 
the white arrow indicates its direction of motion (|Dodson et al. 20031 ). The fact 
that the pulsar is neither at nor moving away from the PWN's center indicates 
that a reverse shock interaction has taken place, (b) T he composite SNR G327.1 - 
1.1. An 843 MHz Molonglo image is shown in red ( Whiteoak fc Green 199fil ). 
while a 0.2-12 keV XMM-Newton image is in blue. The radio morphology con- 
sists of a faint shell enclosing a central PWN. The peak of X-ray emission indicates 
the likely position of an (as yet undetected) pulsar. The offset between the X-ray 
and radio nebulae indicates that the radio nebula is a "relic PWN" as discussed 
in ^13.31 The pulsar is likely to be still moving subsonically through the SNR 
interior, generating a new PWN as it moves away from its birthsite. 
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A4 as 3.1 at this point ( van der Swaluw et al. 2003T ). The injection of energy 
from the pulsar may brighten and re-energize the SNR shell during its p assage 
( Shull. Fesen fc Saken 19891 : Ivan der Swaluw. Achterberg fc Gallant 20021 ). 



3.4 A Pulsar in Interstellar Gas 

Once outside its SNR, a pulsar's motion is often highly supersonic in interstellar 
gas. A bow-shock PWN results, with a potentially large Mach number, Ai S> 1. 

In cases where the pulsar propagates through neutral gas, the forward shock 
driven by the PWN is visible, in the form of Ha emission produced by shock 
excitation and charge exchange (see ^5.2[) . The shocked wind also produces syn- 
chrotron emission, resulting in a bright head and cometary tail, both best seen 
in radio and X-rays (see £ i5.3j) . An example of an interstellar bow shock is the 
structure seen around PSR B 1957+20, show in Figure EJb). 

As the pulsar now moves through the Galaxy, its E drops, and its motion carries 
it away from the denser gas in the Galactic plane where most neutron stars are 
born. Eventually most pulsars will end up with low spin-down luminosities in 





la^VOn* 5S m 30' 56 m OO s Z$ m 3C' 55 m 00' 
Right Ascension (.12000) 



Figure 4: (a) The SNR W44 (G34.7-0.4). T he main panel shows a 1.4 GHz 
VLA image of the SNR (|Giacani et al. 19971 ). while the inset shows 8.4 GHz 
VLA data on th e region surrounding the associated young pulsar B1853+01 
( Frail et al. 199(3 ). whose position is marked by a cross. The pulsar is near- 
ing the edge of the SNR, and now drives a small bow-shock PWN as a result of 
its supersonic moti on, (b) The recycled "Black Widow" pulsar (PSR B1957+20) 
and its bow shock (jStappers et al. 20 03). The green shows Ha emission imaged 
with the Anglo- Australian Observatory, while the red shows X-ray emission ob- 
served with Chandra (the blue emission indicates background stars). The pulsar 
is moving at a position angle of 212° (north through east). 
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low density regions, where they may no longer be moving supersonically and 
their energy output is insufficient to power an observable synchrotron nebula. In 
this final stage, a pulsar is surrounded by a static or slowly expanding cavity of 
relativistic material with a radius > 1 pc and confined by the thermal pressure 
of the ISM ( Arons 19831 . iBlandford et al. 1973h. Deep searches have yet to detect 



such "ghost nebulae" (e.g.. iGaensler et al. 2000l ). 

An alternate evolutionary path may take place for old pulsars in binary sys- 
tems, which can eventually be spun up via accretion from a companion. This 
produces a recycled pulsar with a low value of P but a very rapid spin pe- 
riod, P ~ 1 — 10 ms. Such pulsars can have spin-down luminosities as high as 
E ~ 10 34 — 10 35 ergs s" 1 , which is sufficient to generate observable bow-shock 
nebulae, as shown in Figure Hfb). 

4 YOUNG PULSAR WIND NEBULAE 
4.1 Pulsar Winds 

Despite more than 35 years of work on the formation of pulsar winds, there are 
still large gaps in our understanding. The basic picture is that a charge-filled 
magnetosphere surrounds the pulsar, and that particle acceleration occurs in the 
collapse of charge-separated gaps either near the pulsar polar caps or in outer 
regions that extend to the light cylinder (i.e., to i?LC = c/Vi). The maximum 
potential generated by the rotating magnetic field has been calculated for the 
case of an aligned rotator (i.e., with the magnetic and spin axes co-aligned) by 
Goldreich fc Julian (196^) as: 



2c V10 12 GJ V10 km) Vis/ ' V ; 

where Rns is the neutron star radius. The associated particle current is Nqj = 
(£l 2 B p B? NS ) / Zee where Ze is the ion charge. This current, although considerably 
modified in subsequent models, provides the basis for the pulsar wind. In virtu- 
ally all models, the wind leaving the pulsar magnetosphere is dominated by the 
Poynting flux, F^xB: with a much smaller contribution from the particle energy 
flux, -Fparticie- The magnetization parameter, a, is defined as: 

J= FexB = B 2 

^particle 47T/>7C 2 ' 

where B, p, and 7 are the magnetic field, mass density of particles, and Lorentz 
factor, in the wind, respectively. As the wind flows from the pulsar light cylinder, 
typical values of a > 10 4 are obtained. However, models for t he structure of the 
Crab Nebula ((Kennel fc Coroniti 1984al . lEees fc Gunn 19741 ) require a ^ 0.01 



just behind the termination shock (see ^4.3j) . in order to meet flow and pres- 
sure boundary conditions at the outer edge of the PWN. The high ratio of the 
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synchrotron lumi nosity to the total spin-d own power also requires a particle- 
dominated wind ( Kennel Coroniti 1984bh . and implies 7 ~ 10 6 , a value con - 



siderably higher than that expected in the freely expanding wind ( Arons 2002T ). 
Between the pulsar light cylinder and the position of the wind termination shock 
the nature of the wind must thus chan ge dramatically, although t he mechanism 



tne nature 01 tne wind must tnus mange dramatical ly, altnougn t J 
for this transition is as yet unclear fsee l Arous 20nd .l Melatos 1998) 



The loss of electrons from the polar regions of the star represents a net cur- 
rent that needs to be replenished to maintain charge neutrality This may occur 
through ion outflow in equatorial regions. As discussed in H4A\ ions may con- 
tribute to nonthermal electron acceleration in the inner regions of the PWN. 

4.2 Observed Properties of Young PWNe 

The deceleration of a pulsar-driven wind as it expands into the confines of cold, 
slowly expanding supernova ejecta produces a wind termination shock, at which 
electron/positron pairs are accelerated to ultrarelativistic energies (see El4.3j) . As 
these particles move through the wound-up magnetic field that comprises the 
PWN, they produce synchrotron radiation from radio wavelengths to beyond 
the X-ray band. For a power-law electron spectrum, the constant injection of 
particles p lus a finite synchrotron-emitti ng lifetime lead to a spectral break at a 
frequency ( Ginzburg fe Svrovatskii 19651 ): 



Ub = l0 » ( ^ r v * — r 2 hz, (12) 

V10- 6 Gy V10 3 years/ ' V ; 

where -Bpwn is the nebular magnetic field strength. Particles radiating at fre- 
quencies beyond i/f, reach the outer portions of the PWN in ever-diminishing 
numbers; most radiate their energy before they are able to travel that far. The 
result is that the size of the PWN decreases with increasing frequency, as is clearly 
observed in the Crab Nebula (Fig. For PWNe with low magnetic fields, the 
synchrotron loss times are longer and there may not be a significant difference in 
size between the radio and higher frequency bands (e.g., 3C 58 in Fig. |SJ). 

The morphology of a young PWN is often elongated along the pulsar spin axis 
due to the higher equatorial pressure associa ted with the toroidal magnetic field 
(Bc gelman <fe Li 1991 Ivan der Swaluw 200ah . This effect is seen clearly in many 



PWNe (e.g., Figs. n&El) an d allows one to infer the likely projected orientation 
of the pulsar. As the nebula expands (see ^3.1(1 . Rayleigh- Taylor instabilities form 
as the fast-moving relativistic fluid encounters and accelerates slower-moving un- 
shocked supernova ejecta. These form dense, finger-like filamentary structures 
that suffer photoionization from the surrounding synchrotron emission and ra- 
diate recombinatio n lines in the optical and ultraviolet (UV) bands (Fig. ^b); 
Hester et al. 199*^ ) . The increased density compresses the magnetic field around 



the filaments, causing enhanced synchrotron emission. One thus observes radio 
structures that correspond to the optical/UV filaments. 

At the core of the PWN lies the pulsar itself. As its free-flowing equato- 
rial wind encounters the more slowly-expanding nebula, a termination shock is 
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Figure 5: Images of the PWN 3C 58 (G130. 7+3.1). (a) Radio synchrotron emis - 
sion from the confined wind, with filamentary structure (|B,evnolds fc Aller 19881 ). 
(b) Infrared synchrotron emission with morphology similar to the radio nebula 
(whose outer contour is shown in green), (c) X-ray synchrotron emission (blue), 
thermal emission (red) f rom shock-heated ejecta, and central torus/jet structure, 
shown expanded in (d) ([Slane et al. 20041 ). Images are shown with north up and 
east to the left. The scale for the figures is indicated by the 2 arcmin scale bar 
except for panel (d), where the 20 arcsec scale bar applies. 

formed. Particles accelerated at the shock form a toroidal structure, while some 
of the flow is collimated along the rotation axis, pos sibly contributing to the 
formation of jet-like structures ( Bogovalov et al. 20051 ) . These structures gen- 



erate synchrotron radiation that is observed most readily in the X-ray band 
(Fig . Ud1), although a toroid is also observed in optical images of the Crab Neb- 
ula (|Hester et al. 19951 ). The emission pattern from jets or ring-like structures, 
as well as the larger scale geometry of the PWN, thus provide an indication of 
the pulsar's orientation. As we discuss in H4.41 the emission structures in the 
post-shock and jet regions provide direct insight on particle acceleration, mag- 
netic collimation and the magnetization properties of the winds in PWNe. In 
addition, for pulsars whose proper motion is known, constraints on mechanisms 
for producing this population's high-velocity birth kicks can be derived based on 
the degree of alignment between the pulsar spin ax is and the direction of motion 
(Lai. Chernoff Cordes 200ll: iNg R.omani 20041). 
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4.3 The Wind Termination Shock 

The highly relativistic pulsar wind and its wound-up toroidal magnetic field in- 
flate an expanding bubble whose outer edge is confined by the expanding shell of 
SN ejecta. As the wind is decelerated to match the boundary condition imposed 
by the more slowly-expanding ambient material at the nebula radius, a wind ter- 
mination shock is formed at the radius, R w , at which the ram pressure of the 
wind is balanced by the internal pressure of the P WN (Figure l^b] ) : 



R w = V /^/(4vro;cPpwN), (13) 

where u> is the equivalent filling factor for an isotropic wind, and "PpwN is the total 
pressure in the shocked nebular interior. Upstream of the termination shock, the 
particles do not radiate, but flow relativistically along with the frozen-in magnetic 
field. At the shock, particles are thermalized and reaccelerated, producing syn- 
chrotron emission in the downstream flow. From estimates of the field strength, 
the observed X-ray emission implies Lorentz factors 7 ;> 10 6 in the shock. 

A reasonable pressure estimate can be obtained by inte grating the broad-band 

spectrum of the PWN, using standard synchrotron theory ( Ginzburg fc Svrovatskii 19651 ). 



and assuming equipartition between particles and the magnetic field. Typical 
"Ppwn and E values yield termination shock radii of order 0.1 pc, implying an 
angular size of several arcsec at distanc es of a few kpc . For the Crab Nebula, the 
equipartition field is -Bpwn ~ 300 /uG dTrimble 19821). and the pressu re reaches 
equipartition at a radius of ~ (5 — 20)i?„„ (|Kennel fc Coroniti 1984al ). The as- 
sociated spin-down luminosity, £ = 5 x 10 38 ergs s , yields R w ~ 4 x 10 17 cm, 
consistent with the position of the optical wisps and the radius of the X-ray 
ring seen in Figure HTd). Similar calculations indicate a much weaker field, 
.BpwN ~ 80/iG in 3C 58 (|Green fc Scheuer 19921 ). yielding a termination shock 
radius similar to the Crab, i?„, y 6 x 10 17 cm, given the smaller E of the pulsar 
( Slane. Helfand fc Murray 20021 L This lower field strength is also consistent with 



the fact that the observed size of 3C 58 is similar in the radio and X-ray bands 

(Fig- ED- 

It must be noted that high re solution X-ray observations of 3C 5 8 (plane. Helfand &; Murray 20021 ) , 
G21. 5-0.9 (ICamilo et al. 2005T ). G292.0+1.8 dHughes et al. 200lh and many other 



young PWNe and SNRs do not reveal directly the ring-like emission that is ob- 
served just outside the termination shock in the Crab pulsar, as seen in Fig- 
ure ^d). Rather, the compact emission around the pulsar appears slightly ex- 
tended (see Fig. [5[d] ) , possibly originating from regions similar to the Crab Neb- 
ula's torus, downstream from the termination shock. However, the extent of such 
emission still provides a lower limit on 'PpwN, as well as an indication of the 
pulsar orientation. 

4.4 Formation of Tori, Jets and Wisps 

The geometry implied by the X-ray morphology of the Crab Nebula (see Fig. ^d] ) 
is a tilted torus, with a jet of material that flows along the toroid axis, extending 
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nearly 0.25 pc from the pulsar. A faint counter-jet accompanies the structure, 
and the X-ray emission is significantly enhanced along one edge of the torus. Both 
effects are presumably the result of Doppler beaming of the outflowing materia l, 



whereby the X-ray intensity, I, varies with viewing angle as ()Pelling et al. 19871 ): 

r+i 



I 

io 



1 — (3 cos (f> 



(14) 



where f3c is the flow speed immediately downstream of the termination shock, (f> 
is the angle of the flow to the line-of-sight, and Io is the unbeamed intensity. 

Similar geometric structures are observed in G54. 1+0.3, for which Chandra 
observations reveal a central point-like source surrou nded by an X-ra y ring whose 
geometry suggests an inclination angle of about 45° ( Lu et al. 2002F ). The X-ray 



emission is brightest along the eastern limb of the ring. If interpreted as the result 
of Doppler boosting, this implies a post-shock flow velocity of ~ 0.4c. The ring 
is accompanied by faint bipolar elongations aligned with the projected axis of 
the ring, consistent with the notion that these are jets along the pulsar rotation 
axis. The total luminosity of these structures is similar to that of the central ring. 
This is to be contrasted with the Crab and 3C 58 (Figs. [l]&0~J), for which the 
torus outshines the jets by a large factor. Moreover, for G54. 1+0.3 the brighter 
portion of the outflow lies on the same side of the pulsar as the brightest portion 
of the ring, which is inconsistent with Doppler boosting. A similarly troubling 
observation is that the brightness distribution around the inner ring of the Crab 
also fails to show Doppler brightening consistent with that seen in its surrounding 
torus; the brightness is reasonably uniform except for some compact structures 
that vary in position and brightness with time (see H4.5|) . 

The formation of these jet /torus structures can be understood as follows. Out- 
side the pulsar magnetosphere, the particle flow is radial. The rotation of the 
pulsar forms an expanding toroidal magnetic field for which the Poynting flux 
varies as sin 2 ip, where ijj is the angle from the rotation axis. Conservation of 
energy flux along flow lines leads to a latitude dependence of the Lorentz factor 
of the wind of the form 7 = 70 + 7m sin 2 -0 ( Boeovalov &: Khangoulvan 20021 ) 



where 70 is the wind Lorentz factor jus t outside the light cylinder (70 ~ 10 2 in 
standard models for pulsar winds; e.g., Daughertv k, Harding 198*21 ) . and 7 m is 



the maximum Lore ntz factor of the preshock wind particles (7 m ~ 10 6 near the 
termination shock; Kennel fc Coroniti 1984ah . From Equation (jllj) . we see that 



this corresponds to a latitude variation in the magnetization parameter also, with 
a much larger at the equator than at the poles. This anisotropy results in the 
toroidal structure of the downstream wind. Moreover, modeling of the flow con- 
ditions across the shock shows that magnetic collimation produces jet-like flow s 
along the rotation axis ( Bogovalov et al. 20051 . Komissarov fc Lvubarskv 20041 ). 



This collimation is highly dependent on the magnetization of the wind. For 
a ^ 0.01, magnetic hoop stresses are sufficient to divert the toroidal flow back 
toward the pulsar spin axis, collimating and a ccelerating the flow to speeds of 
~ 0.5c ( Del Zanna. Amato fc Bucciantini 20041 ): smaller values of a lead to an 
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increase in the radius at which the flow is diverted. The launching points of the 
jets observed in the Crab Nebula, which appear to form much closer to the pulsar 
than the observed equatorial ring, apparently reflect such variations in the value 
of a. Near the poles a is large, resulting in a small termination shock radius and 
strong collimation, while near the equator a is much smaller an d the termination 
shock extends to larger radii ( Bogovalov fc Khangoulvan 200 j ). 



Because 7 is smaller at high latitudes, models for the brightness of the jet-like 
flows produced by the collimation process fall short of what is observed. However, 
kink instabilities in the toroidal field may transform magn etic energy into particl e 
energy, thus accelerating particles and brightening the iet feogovalov et al. 20051 ) . 



Such instabilities limit the duration of the collimation (|Begelman 19981 ) . Evidence 
for the effects of kink instabilities is seen in the curved nature of many PWN jets, 
particularly for the Vela PWN where the iet morphology is observed to change 



dramatically on timescales of months ((Pavlov et al. 200.1l ). There appears to be 



a wide variation in the fraction of E channeled into PWN jets, ranging from 
~ 2.5 x 10" 5 in 3C 58 to nearly 10~ 3 for the PWN powered by PSR B1509 - 
58, based on their synchrotron spectra ( Gaensler et al. 2002al . Slane et al. 20041 ). 



This apparently indicates considerable differences in the efficiency with which 
additional acceleration occurs along the jets. 

The torus surrounding the Crab pulsar is characterized by the presence of 
wisp-like structures (see Fig. [Jd] ) , whose position and brightness vary with time 
in the optical, infrared and X-ray bands, and which emanate from the termina- 
tion sh ock and move outward through the tor us with inferred outflow speeds of 
0.5c dHester et al. 200l iMelatos et al. 2005h. The exact nature of these struc- 



tures is not fully understood. iHester et al. f 20021 ) suggest that they are formed 
by synchrotron instabilities. However, the position of the arc-like structure sur- 
rounding PSR B1509-58 seems inco nsistent with this hyp othesis given the much 
lower magnetic field in the PWN (jGaensler et al. 2002aT l . An alternative sug- 



gestion is that the wisps are the sites of compression of the electron/positron 
pair plasma on scales of the cyclotron gyration radius of ions in the outflow, 
which can also explain the radius of the X-ray ar c seen around PSR B 1509-58 
( Gallant &: Arons 19941 . Spitkovskv &: Arons 2004, and references therein). 



It is worth noting that VLA observations of the Crab Nebula show variable 
radio structures very similar to the optical and X-ray wisps, indicating that ac- 
celeration of the associated par ticles must be occurrin g in the same region as for 
the X-ray-emitting population (jBietenholz et al. 20041V 



4.5 Filamentary and Compact Structures in PWNe 

In the Crab Nebula, an extensive network of filaments is observed in Ha, [O ill] 
and other optical lines, surrounding the nonthermal optical emission from the 
PWN (Fig.^b]). The detailed morphology and ionization structure of these fila- 
ments indicate that they form from Rayleigh- Taylor instabilities a s the expanding 
relativistic bubble sweeps up and accelerates slower moving ejecta ( Hester et al. 199fil ). 
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Magnetohydro dynamic (MHD) simulations support this picture, in dicating that 

60%-75% of the swept-up mass can be concentrated in sucn filaments feucciantirn et al. 2004 
Jun 1998ft . As the expanding PWN encounters these filaments, compression in- 



creases the density and magnetic field strength, forming sheaths of en hanced syn- 
chrotron emission observed as a corresponding shell of radio fila ments dE,evnolds 1988ft . 
X-ray observations reveal no such filaments in the Crab Nebula ((Weisskopf et al. 20001 ) . 
This is presumably because the higher energy electrons required to produce the 
X-ray emission suffer synchrotron losses before reaching the outer regions of the 
PWN, consistent with the smaller extent of the nebula in X-rays than in the 
radio. 

A different picture is presented by X-ray observations of 3C 58, which reveal 
a complex of loop-like filaments most prominent near the central region s of the 
PWN, but evident throughout the nebula (Tig. IBTc): ISlane et al. 2004ft . These 
nonthermal X-ray struc tures align with filaments observed in the radio band 
( Reynolds Alter 1988ft . Optical observations of 3C 58 also reveal faint filaments 
( van den Bergh 19781 ) . whose origin is presumably similar to those in the Crab 
Nebula. However, a detailed X-ray/optical comparison shows that most of the 
X-ray filaments do not have optical counterparts. While comparisons with deeper 
optical images are clearly needed, the fact that many of the X-ray features without 
optical counterparts are brighter t han average sugges ts that these may actually 
arise from a different mechanism. Slane et al. r2004ft propose that the bulk of 
the discrete structures seen in the X-ray and radio images of 3C 58 are magnetic 
loops torn from the toroidal field by kink instabilities. In the inner nebula, the 
loop sizes ar e similar to the s ize of the termination shock radius (~ 0.1 pc), as 
suggested by Begelman f 1998ft . As the structures expand, they enlarge slightly 
as a consequence of decreasing pressure. 

There is also considerable loop-like filamentary structure evident in hi gh reso- 
lution X-ray images of the Crab Nebula fFig. OTd^: IWeisskopf et al. 2000ft . These 
filaments appear to wrap around the torus, perpendicular to the toroidal plane, 
and may be signatures of kink instabilities in the termination shock region. 

In some PWNe, compact knot-like structures are observed close to the pulsar , 
which dissipate and reappear on timescales of order months ( Hester et al. 200^ . 
Melatos et al. 2005ft . Examples can be seen for t he Crab in Figure HTd) , and sim- 



ilar structures are also seen for PSR B1509-58 (iGaensler et al. 2002aft . Some of 
these features appear in projection inside the termination shock region. However, 
it is believed that they actually correspond to unstable, quasi-stationary shocks 
in the region just outside the termination shoc k, at high latitudes where the sho ck 
radius is small due to larger values of a (e.g., Komissarov fc Lvubarskv 20041 ). 



4.6 PWN Spectra 

As noted in ^2.21 the spectra of PWNe are characterized by a flat power law 
index at radio wavelengths {a ~ —0.3) and a considerably steeper index in X- 
rays (T k, 2). The nature of this spectral steepening is not understood. Simple 
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assumptions of a power-law particle spectrum injected by the pulsar would predict 
a power-law synchrotron spectrum, with a break associated with the aging of the 
particles (see Eqn. [T2] the expected increase in spectral index is Aa = 0.5, 
which is smaller than what is typically observed ( Woltier et al. 199?ft . Moreover, 



for many PWNe a change in spectral index is inferred at low frequencies tha t 
would imply unrealistically high magnetic fields (e.g., Green &: Scheuer 1992F ). 



Relic breaks in the spectrum can be produced by a rapid decline in the pulsar 
outp ut over time, and these breaks propagate to lower frequencies as the PWN 
ages ( Pacini fc Salvati 19731 ). The inherent spectrum of the injected particles, 



which may deviate from a simple power law, as well as modifications from discrete 
acceleration sites, all contribute to a complicated integrated spectrum. As a 
result, the interpretation of spectral steepening as being due to synchrotron losses 
can lead to drastically wrong conclusions about PWN properties. 

At frequencies for which the synchrotron lifetime is shorter than the flow time 
to the edge of the PWN, one expects a steepening of the spectrum with ra- 
dius. Radial steepening is ind eed observed in the X-ray spectra of the Cra b 
Nebula, 3C 58 and G21.5-0.9 (jSlane et al. 200(1 120041 IWillingale et al. 200lh . 



but the spectra steepen rather uniformly with radius whereas generalizations of 
the lKennel k. Cororiiti ri984tJl model predict a much more rapid steepening near 
the outer regions (|Revnolds 20031 ). Some mixing of electrons of different ages at 
each radius seems to be required, perhaps due to diffusion processes in the PWN. 

5 BOW SHOCKS AROUND HIGH VELOCITY PULSARS 

Pulsars are born with high space velocities, typically Vpsr = 400-500 km s _1 , 
but for some sources exceeding 1000 km s" 1 . These high velocities are almost 
certai nly the re sult of kicks given to the star during or shortly after core col- 
lapse ( Lai 200-ih . Young pulsars thus have the highest velocities of any stellar 



population, and many have sufficient speeds to eventually escape the Galaxy. 

As discussed in ^13.31 a pulsar's ballistic motion allows it to eventually escape its 
original PWN, and to propagate through the shocked ejecta in the SNR interior. 
At first the pulsar's motion will be subsonic in this hot gas, but by the time 
the pulsar nears the edges of the SNR, the sound speed drops sufficiently for the 
pulsar's motion to be supersonic. In the simplest situation of a spherical SNR 
in the Sedov phase, expanding into a uniform medium, this transition occurs at 
half the crossing time (given in Eqn. ^5] below), at which point the pulsar has 
traveled 68% of the distance from the center of the SNR to its edge (van der 
Swaluw et al. 2004). This simple result is independent of Vpsr, no or E$n- The 
pulsar's supersonic motion now produces a PWN with a bow shock morphology. 

Since the SNR is decelerating, the pulsar ultimately penetrates and then es- 
capes the shell. A pulsar moving at Vpsr > 650 km s -1 will escape wh ile the 
SNR is still in the Sedov phase, after a time ( van der Swaluw et al. 20031 ): 



/ E SN \ 1/3 ( n X- 1 / 3 / Vpsr \ ~ 5/3 



t cross = 44 —^3 — coa kyr. (15) 

V10 5i ergs/ VI cm 6 J V500 km s 
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At times t > t cross , a pulsar proceeds to move through the ambient ISM. 

The speed of sound in interstellar gas is a a function of temperature: typical 
values are approximately 1, 10 and 100 km s _1 for the cold, warm and hot 
components of the ISM, respectively. Thus except in the case of a particularly 
slow moving pulsar moving through coronal gas, a pulsar will move supersonically 
and drive a bow shock through the ISM. 



5.1 Theoretical Expectations 

The pulsar wind in a bow shock is decelerated at a termination shock, just as for 
younger PWNe (see t !4.3J) . However, the external source of pressure balance is now 
ram pressure from the pulsar's motion, rather than the internal pressure of the 
shocked wind. Furthermore, since ram pressure is not isotropic, the termination 
shock radius varies as a function of angle with respect to the pulsar's velocity 
vector. In the direction of the star's motion, the termination shock radius is 
referred to as the "stand-off distance", R w o, and is defined by (cf. Eqn. [13]): 

Po^sr, (16) 



where po is the ambient density. If the wind is isotropic and u = 1, then at a polar 
angle 6 with respect to the bow shock's symmetry axis, t he analytic so lution for 
the termination shock radius as a function of position is ( Wilkin 199fil ): 



R w (9) = R w0 csc 9 ^3(1-6 cot 6). (17) 

It is important to note that the above solution assumes an efficiently cooled 
thin-layer shock, in contrast to the double shock expected for pulsar bow shocks. 
Full hydrodynamic and MHD simulations show that Equation (|17|) is a reasonable 
approximatio n in regions near the apex (6 <; lp, but performs more poorly further 
downstream ( Bucciantini 200211 van der Swaluw et al. 20031 ) . 



A result of such simulations is shown in Figure a). The double-shock struc- 
ture is clearly apparent, consisting of a forward shock where the ISM is heated, 
plus the termination shock where the pulsar's wind decelerates. As expected, the 
termination shock is not of uniform radius around the pulsar: specifically, for low 
Mach numbers, A4 ~ 1 — 3 (as may be appropriate for pulsars traveling super- 
sonically inside their SNRs; see H3.3l and van der Swaluw et al. 2004), the ratio of 
ter mination shock radi i between polar angles 6 = t t and 8 = is approximately 
M (|Bucciantini 2002al . Ivan der Swaluw et al. 20031 ) , but for M > 1 (typical of 



bow shocks in the ambient ISM; see ^3.4J) , this ratio approach es a limit of ~ 5 — 6 
( Bucciantini. Amato fc Del Zanna 200,4 iGaensler et al. 200 J ). 



5.2 Observations: Forward Shock 

If a pulsar drives a bow shock through neutral gas, then collisional excitation and 
charge exchange occur at the forward shock, generating optical emission in the 





Figure 6: (a) A hydro-dynamic simulation of a pulsar bow shock, adapted from 
Gaensler et al. (2004). The pulsar, whose position is marked with an asterisk, 
is moving from right to left with a Mach number A4 = 60. The intensity in 
the image and the scale-bar indicate density, in units of log 10 (/9o/10 -24 g cm -3 ), 
(b) Chandra X-ray (blue) and VLA radio (red) images of G359.23-0.82 ("th e 
Mouse"), the bow shock associated with PSR J1747-2958 (|Gaensler et al. 20041 ). 
The white arrow marks a bright compact region of X-ray emission behind the 
apex, which possibly corresponds to the surface of the termination shock. 



Balmer lines ( Bucciantini 2002bl . Bucciantini fc Bandiera 200ll ). Indeed several 



pulsar bow shocks have been identified in th e 656-nm Ha line, associa ted with (in 
order of disco very): B1957+20 (Fig. ETV): iKiilkarni Hester 1988ft. B2224+6 5 
( "the Oi]itar": ICordes. Romani fc Lnndgren 19931) . .70437-471 5 dBell et al. 199,fll. 



RX .T1 856.5-3754 (Ivan Kcrkwiik .V Kiilkarni 200lh. B 0740-28 d.lones. Stunners .V Gaensler 2002h 



and J2124-3358 (|Gaensler. Jones Stanners 2002bh . 



If the distance to the system is known, R w q can be directly measured, provided 
that one adopts a scaling factor of ~ 0.4 — 0.6 to translate be tween the observed - 
radius of the forward shock to that of the termination shock ( Bucciantini 2002"al . 



van der Swaluw et al. 200*31 ). If E and Vpsr have been measured, Equation (fl"6j) 



can then be applied to yield p$. This is an approximation, since u = 1 is usually 
assumed, and because of the unknown inclination of the pulsar's motion to the 
line of sight, 3 but certainly suggests ambient number densities ~ 0.1 cm" 3 , as 
expected for warm neutral ambient gas (Chatterjee & Cordes 2002, Gaensler et al. 
2002b). For pulsars for which Vpsr is not known, one can write: 

PoVpsr = Ti-A^ism, (18) 

where 71 = 5/3 and V\sm are the adiabatic coefficient and pressure of the ISM, 
respectively. Estimates for T^ism and to in Equations ((TB|) & (fTH)) then yield A4. 

For the bow shocks associated with RX J1856.5-3754, PSR J0437-4715 and 
PSR B1957+20 (Fig. Bib]), the shape of the forward shock is a good match to 

3 Note that the correction factor for motion inclined to the line of sight cannot be derived 
through simple trigonometry (Gaensler et al. 2002b). 
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the solution predicted by Equation (|17|) and simulated in Figure a). However, 
the optical emission round PSRs B2224+65, J2124-3358 and B0740-28 all show 
strong deviations from the expected shape, in that there are abrupt kinks and 
inflection points in the Ha profile (Jones et al. 2002, Chatterjee & Cordes 2002). 
Furthermore, in the case of PSR J2 124-3358 there is an apparent rotational 
offset between the symmetry axis of the bow shock and the velocity vector of 
the pulsar (Gaensler et al. 2002b). These systems imply the presence of some 
combination of anisotropies in the pulsar wind (as are observed in young pulsars; 
see t )4,4[) . gradients and fluctuations in the density of the ISM, or a bulk velocity 
of ambient gas with respect to the pulsar's local standard of re st. Extensions o f 
Equatio n (1171) to accou nt for these effects have been presented by Bandiera (19931 ) 
and by IWilkin (2000l ). and have been applied to interpret the morphology of 
PSR J2124-3358 by Gaensler et al. (2002b). 

5.3 Observations: Termination Shock 

Just as for the PWNe around the youngest pulsars discussed in M.'Al particles in 
the pulsar wind inside a bow shock will be accelerated at the termination shock, 
producing non-thermal synchrotron emission that can be potentially observed in 
the radio and X-ray bands. Indeed cometary radio and X-ray PWNe aligned 
with the direction of motion has been now identified around m any pulsars, con- 
vincing examples of which in clude PSRs B1853+01 ( Fig. BTa): iFrail et al. 199fil : 
Petre. Kuntz Shelton 2002h , B1957+20 (Fie . HTM: 1st minors <>l nl. 2003!)lmd 



B1 757-24 ("the Duck" : IFrail fe Kulkarni 1 99 ifTknsni et «l. 200 lh. 



The most spectacular example of this class is G35 9. 23-0.82, the X-ray /radio 

bow shock powered by PSR J1747-2958 ( "the Mouse" : iGa^n si er et al. 20ollY 1 isef-Zadeh fe Rally 1987ft 



multi- wavelength observations of which are shown in Figure Efb). The extent of 
the radio trail in this system is larger than in X-rays, reflecting the difference in 
synchrotron lifetimes between these bands (cf. H4.2JI . 

The X-ray morphology of the Mouse in Figure Efb) appears to consist of two 
components: a bright compact region extending ~ 0.2 pc from the pulsar, su- 
perimposed on a larger fainter component extending ~ 1 pc from the pulsar. 
Comparison with the hydrodynamic simulation in Figure a) suggests that the 
bright component of the X-ray emission corresponds to the surface of the wind 
termination shock, while fainter, more extended X-ray emission originates from 
the shocked wind ( Gaensler et al. 2004ft . In this identification, the bright X-ray 



component of this bow shock corresponds to the inner toroidal rings discussed in 
^14.41 but elongated due to the pulsar motion. This interpretation is supported 
by the fact that the extent of the termination shock region along the symmetry 
axis (in units of R w q) should be a function of Mach number, as w as discussed 
in M5.ll Indeed for the Mouse, with M ~ 60 ((Gaensler et al. 20o3 ) . this bright 
region is about twice as long (relative to -R „,n) as for the pulsar bow shock seen 



inside the SNR IC 443 (jOlbert et al. 200 lh . whic h must have M > 3 since it is 



moving through shocked gas in the SNR interior ( van der Swaluw et al. 20031 ) . 
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In PWNe for PSR s B1757-24 flKaspi et al. 200lft and B1957+20 (Fig. Efb) 
Stappers et al. 2003ft . only a short (~ 0.1 pc) narrow X-ray trail is apparent. Such 



features have been int erpreted as a rap i d back-flow or nozzle, which transports 
particles downstream ([Kaspi et al. 200ll : IWang. Li fc Begelman 1993ft . However, 



Figure ©suggests that the short trail seen behind PSRs B1757-24 and B1957+20 
is the surface of the termination shock, a nd that emission from the post-shock 
wind further downstream is too faint to see ( Gaensler et al. 20041 . Gvaramadze 2001 ) . 



Deeper X-ray observations are required to test this possibility. 

Hydrodynamic models predict that the pulsar's motion should divide the post- 
shock emitting region of a bow shock into two distinct zones: a highly magnetized 
broad tail originating from material shocked at 9 <^ ^, plus a more weakly magne- 
tized, narrow, collimated tail, produced by material flowin g; along the axis 9 ~ tc 
( Bucciantini 2002al : Romanova. Chulskv fc Lovelace 200$ ) . These two structures 



are both apparent in the Mouse (|Gaensler et al. 20041 ). Through relativistic MHD 
simulation, this and other issues related to the structure of the post-shock flow 
are now being explored (e.g., Bucciantini et al. 2005). 

6 OTHER TOPICS AND RECENT RESULTS 
6.1 Pulsars and PWNe in Very Young SNRs 

The youngest known Galactic PWNe are the Crab Nebula and 3C 58, powered by 
pulsars thought to correspond to the SNe of 1054 CE and 1181 CE, respectively. 4 
It would be of great interest to identify PWNe at earlier evolutionary stages. 
SN 1987A formed a neutron star, but deep searches ha ve failed to detect a 



central object, down to a luminosity <^ 10 34 ergs s _i (e.g.. iGraves et al. 
This is well below the luminosity of the Crab Nebula, and may indicate that the 
central neutron star has collapsed further int o a black hole, accretes from f all-back 



material, or does not generate a wind (e.g., Frver. Colgate &; Pinto 19991 ) . 



Searches for PWNe in other extr agalactic SNe and SNRs have generally not 
produced any convincing candidates dBartel fc Bietenholz 20051 . iRevnolds fc Fix 1987ft . 



However, recent high resolution radio images have revealed the gradual turn-on 
of a central flat-spectrum radio nebula in SN 1986J, which may co rrespond to 
emission from a very young PWN dBietenholz. Bartel fc Rupen 20041 ). New wide- 



field radio and X-ray images of other galaxies may lead to further identification 
of young PWNe, while optical and UV spectroscopy of recent SNe may identify 
PWNe through emissio n lines which broaden with t ime as gas is swept up by the 
expanding pulsar wind ( Chevalier &: Fransson 1992ft . 



4 The association of 3C 58 with 118 1 CE is not completely secure; see 
IStephenson fc Green (2002T) and lChevalier (20051) . and references therein. 
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6.2 Winds from Highly Magnetized Neutron Stars 

A growing population of neutron stars have surface magnetic fields (inferred via 
Eqn. 0]) above the quantum critical limit of B p = 4.4 x 10 13 G. The proper- 
ties of these stars indicate that they ar e comprised of two app arently distinct 
populations: the high-field radio pul sars dMcLaughlin et al. 20031 . and references 
therein), and the exotic magnetars ( Woods fc Thompson 2006f ). The winds and 



PWNe of these sources potentially provide a view of different spin-down processes 
than those seen in normal pulsars. 

Most of the high-field radio pulsars have E ~ 10 32 — 10 34 ergs s . For typ- 
ical efficiency factors r\R,r\x <^ 10~ 3 , this implies PWNe too faint to be de- 
tectable. However, the very young pulsars J1846-0258 (B p = 4.9 x 10 13 G) 
and J1119-6127 (B p = 4.1 x 10 14 G) have high spin-down luminosities (E > 
10 36 ergs s" 1 ) and are near the centers of SNRs. In both cases PWNe are de- 
tected, although with very different properties. PSR J1846-0258 puts a large 
fraction {r]x ~ 0-2) of its spin-down power int o a luminous X-ray PWN ~ 
2 pc in extent (|Helfand. Collins fc Gotthelf 200.^ 1 . while PSR J1119-6127 pow- 



ers an under-luminous (vx ~ 2 x 10 4 ) and small (~ 0.2 pc) X-ray nebula 
( Gonzalez fc Safi-Harb 20081 ). Clearly PWN properties are dominated by fac- 



tors such as age, environment and evolutionary state (see ^SJ rather than the 
associated pulsar's surface magnetic field. 

Magnetars also spin down, albeit in some cases not smoothly ( Woods et al. 20021 ). 



Just as for radio pulsars, this rotational energy output is thought to go into a rel- 
ativistic wind, but traditional PWNe have not been detected around magnetars, 5 
presumably because E £ 10 34 ergs s _1 for all these sources. Magnetars likely 
experience an enhanced torque over the dipole spin-down presumed to act in ra- 
dio pulsars, as a result of either Alfven waves and outflowing relativistic particles 
driven by seismic activity, or by a large-scale twist of the external magnetic field 
( Harding. Contopoulos fc Kazanas 19991 : Thompson. Lvutikov fc Kulkarni 2002h . 



Under either circumstance, t he spin-down behavior deviates from that described 
in Equations © to (e.g.. iThompson et al. 2000h . 



A transient radio PWN was proposed to account for the short-lived radio 
nebula seen in 1998 following t he giant flare from the magnetar SGR 1900+14 
( Frail. Kulkarni fc Bloom 1999T l. However, recent observations of a radio nebula 



in the aftermath of a giant flare from SGR 1806-20 suggest that the synchrotron 
emission from these nebulae are powered by ejected bar yonic material, so th at 
these sources are more analogous to SNRs than to PWNe (|Gaensler et al. 2005^1 . 6 



5 The magnetar SGR 1806-20 was originally presumed to power the rad io nebula G10.0-0.3 , 
but a revision in the position of this neutron star now makes this unlikely (iHurlev et al. 19991) . 

6 A possibly transient radio source has also recently been identified coin cident with the flarin g 
magnetar XTE J1810-197, but the nature of this source is as yet unclear ijHalnern et al. 20051) . 
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Rifjht Ascension (J2OO0) 



Figure 7: Multi-wavelength images of the PWN powered by the young pulsar 
B1509-58. ROSATPSPC data (in blue contours, at levels of 5%, 10%, 20%, 40% 
and 60% of the peak) show the extent of the X-ray PWN ( Trussoni et al. 19961 ) , 
while 843 MHz Molonglo data (in r ed) correspond to the surrounding SNR 
G320. 4-1.2 (IWhiteoak k Green 199fih . TeV emission from HESS is shown in 
green ( Aharonian et al. 2005bl ). The cross marks the position of PSR B1509-58 



6.3 TeV Observations of PWNe 



The Crab Nebula is a well-known source of TeV gamma-rays dWeekes et al. 19891 . 
and references therein). This emission is well explained as inverse Compton 
(IC) emission, the relativistic particles in the shocked wind acting as scattering 
cente rs for the synchrotron pho tons that they themselves emit at lower ener- 
gies (|Atovan k Aharonian 199~6h . Under this interpretation, the emitted spec- 
trum can be modeled to provide the mean and spatial distribution of the neb- 
ular magnetic field strength, and hence the PWN's particle content, the time- 
averaged injection rate of part icles, and an indepe ndent estimate of the mag- 



netization parameter, a (e.g., de Jaeer et al. 1996h . The estimated values of 



-B pwn and a are in good ag reement with those derived from the MHD model 
of lKennel k Coroniti Cl984ah . 

The new generation of ground-based Cerenkov detectors (most notably the 
High Energy Stereoscopic System, HESS) have now begun to detect other PWNe 
in the TeV band. These detections indicate that acceleration of particles to 
considerable energies must have occurred, and provide estimates of the nebular 
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magnetic field strength which can be used in modeling and interpreting the other 
nebular structures discussed in fjl| The much lower synchrotron luminosities of 
these other sources compared to the Crab Nebula implies that the seed photons 
for IC scattering are in these cases primarily external, originating from a combi- 
nation of the cosmic microwave background and a local contribution from dust 
and starlight. Recent TeV de tections of PWNe by H ESS include GO.9+0.1 and 



G320.4-1.2 / PSR B1509-58 (|Aharonian et al. 2005al lhh. As shown in Figured 



the latter is spatially resolved by HESS and has a TeV morphology which is a 
good match to the X-ray synchrotron nebula. Such observations can potentially 
provide direct measurements of spatial variations in the magnetic fields of PWNe. 

6.4 Pulsar Winds in Binary Systems 

The recently discovered dual-line double pulsar PSR J0737-3039 consists of a 
23-ms p ulsar ( "A" ) and a 2.8-s pulsar ("B") in a 2.4-hour orbit, viewed virtually 
edge-on ( Lvne et al. 2004| ) . This system is proving to be a remarkable new probe 



of pulsars and their winds, providing information at much closer separations to 
the pulsar than is possible for the sources discussed in £@] & 

The line of sight to pulsar A passes within ~ 0.01 light-seconds of pulsar B, 
well within the unperturbed light cylinder radius of the slower pulsar. For ~ 30 s 
at conjunction, the pulse-averaged radio emission from A is modulated with a 
complicated time-dependence, showing inter mittent periodicities at bo th 50% and 
100% of the rotational period of pulsar B ( McLaughlin et al. 2004bl ). Detailed 



modeling shows that this behavior can be interpreted as synchrotron absorption 
from a relativistic plasma confined by a dipolar magnetic field, providing direct 
evidence for the field geome t ry commonly adopted in pulsar electrodynamics 
( Lvutikov fc Thompson 200,1 iRafikov fc Goldreich 2005h . 



Drifting sub-pulses in the pulsed emission from pulsar B are also observed, 
with fluctuations at the beat frequency between the periods of A and B, and 
with a separation betwee n drifting features corresponding to the period of A 
( McLaughlin et al. 2004^ . This provides clear evidence that pulsed radiation 



from A is interacting with the magnetosphere of B, supporting the picture that 
pulsar winds are magnetically dominated (a 1) in their inner regions, as dis- 
cussed in ij4.ll 

Other binary systems also provide information on conditions very close to the 
pulsar. PSR B1957+20 is in a circular 9.2-hour orbit around a ~ 0.025 M Q com- 
panion, interactio n with which produc es a termination shock just 1.5 x 10 11 cm 
from the pulsar 7 ( Phinnev et al. 19881 ). The X-ray flux from t his nebula shows 



possible modulation at the orbital period ( Stappers et al. 2003h. as should resul t 



from Doppler boosting of the flow around the companion (lArons fc Tavani 1993h . 
PSR B1259-63 is in a highly eccentric 3.4-year orbit around a Be star. Near pe- 
riastron, the pulsar is subject to a time-varying external pressure, producing a 
transient X-ray/radio synchrotron nebula plus TeV emission from IC scattering of 

7 On much larger scales, PSR B1957+20 also powers a bow shock, as shown in Figure |ljb). 
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light from the companion star (e.g., Aharonian et al. 2005cl . Johnston et al. 20051 ). 
At future periastra of this system, coordinated X- and 7-ray observations with 
INTEGRAL, GLAST and HESS can directly probe particle acceler ation in this 



pulsar's wind ((Kirk. Ball Sz Skiaeraasen 19991 : iTavani &: Arons 19971 ) 



SUMMARY POINTS 

1. A magnetized relativistic wind is the main reservoir for a pulsar's rotational 
energy loss. The termination of this wind due to surrounding pressure 
produces a pulsar wind nebula, usually observed as a synchrotron nebula. 
A high velocity pulsar can also produce a line-emitting optical bow shock 
where the pulsar wind shocks surrounding gas. 

2. Pulsar wind nebulae move through a series of distinct evolutionary states, 
moderated by the pulsar's location (inside a SNR vs in interstellar gas), the 
ambient conditions (cold ejecta vs shocked ejecta vs ISM) and the Mach 
number of the pulsar (subsonic vs mildly supersonic vs highly supersonic). 

3. High resolution X-ray observations of young pulsars reveal the imprint of 
the rotation axis on the morphology of the surrounding PWN, in the form 
of equatorial tori, polar jets, and overall elongation of the nebula. Using 
these structures, one can locate the wind termination shock and can infer 
the composition, flow speed and geometry of the pulsar's wind. 

4. An increasing number of bow shocks are being found around high-velocity 
pulsars. These systems impose a second axis of symmetry on the PWN, 
providing additional probes of the pulsar's wind and environment. 
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